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I. INTRODUCTION 

The purpose of this study is to examine the effectiveness of the 
geometrical theory of diffraction (GTD) in calculating the radiation 
from slots in perfectly-conducting surfaces which are partially 
covered with dielectric. Configurations of this type are of interest 
in the design of antenna systems for spacecraft. We consider the 
radiation from a narrow slot in a dielectric-covered perfectly- 
conducting surface terminated at an edge as shown in Fig. 1. This is 
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P_: The Field Point in the 

Shadow Region 


S^: Geometrical Optics Ray 

S-| : Singly-diffracted Ray 

S2^ Doubly-diffracted Ray 

S Surface Ray Traveling 
Parallel to the 
Dielectric Sheet 


Fig. 1 



one of the most elementary structures involving dielectric-covered 
slots; furthermore, measured patterns are available to check the GTD 
analysis [l]. 

A theoretical approach to the problem based on the geometrical 
theory of diffraction [2] is described. The effect of the truncation of 
the dielectric cover has been ignored in previous analysis, yet this 
yields the diffracted ray which plays a dominant role in determining 
the radiation pattern. The total far-zohe field is composed of a 
geometrical optics field and a diffracted field. The geometrical 
optics field is the direct radiation from the slot to the field point 
P; it vanishes in the shadow region below the plane containing the 
surface AB. The slot also excites surface waves which are incident 
at the termination of the dielectric cover at A and B. The termi- 
nations yield singly-diffracted rays and reflected surface waves. 

For simplicity the thickness t of the dielectric cover is restricted 
so that only the dominant TMq surface wave mode is excited. The 
contributions from the geometrical optics ray, the singly- diffracted 
rays and all significant multiply-diffracted rays are summed to give 
the field at P. ' . 

The geometrical optics field and the field of the incident 
surface wave are obtained from the solution to the problem of a narrow 
slot radiating through a dielectric-covered ground plane of infinite 
extent. The diffraction coefficient and reflection coefficient for 
the surface wave incident on the edge A (or B) are determined in the 
following way. The surface impedance of the grounded, dielectric cover 

kp 

(1) Zs ' zo r 
0 

where k^ and Zq are the wave number and characteristic impedance of 
free space, and k^, which is found from thesolution of a transcen- 
dental equation, is the wave number of the surface wave in the direction 
normal to the dielectric-covered surface (kp is imaginary for a lossless 
dielectric cover). Next, consider a right-angle wedge, one of whose 
faces is perfectly-conducting, the other with surface impedance Zg. The 
canonical problem of the diffraction of a surface wave by this right- 
angle wedge [3] is used to find diffraction and reflection coefficients, 
which are adequate for sufficiently thin dielectric covers. However, 
as the thickness of the cover increases, the radiation from the 
vertical end face AA' (or BB') of the dielectric cover cannot be 
neglected. This is taken into account using a Ki rchhoff-type approxi- 
mation, where the total electric field within the dielectric cover at 
its end face radiates in the presence of a right-angle wedge. This 
contributes a second term to the diffraction coefficient previously 
obtained. 
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The calculated and measured patterns compare well, which justifies 
the use of the geometrical theory of diffraction in this type of 
problem together with the approximations which have been used to 
find the diffraction coefficient. The patterns are very sensitive 
functions of the parameters, such as the thickness of the dielectric 
cover and the frequency. When the dielectric cover is lossy, the 
contributions from the diffracted rays decrease as the loss tangent 
of the dielectric cover increases. This is due to the fact that the 
surface wave is attenuated as it is incident at the termination of 
the dielectric cover. The ripples in the patterns in the illuminated 
region diminish and the level of radiation in the ' shadow region 
is reduced. 

II. METHOD OF ANALYSIS 

According to the geometrical theory of diffraction, the total 
far-zone field is composed of a geometrical optics field and a dif- 
fracted field. In the problem of a narrow slot in a dielectric- 
covered perfectly-conducting surface terminated at an edge (see 
Fig. 1), the geometrical optics field is the direct radiation from 
the slot to the field point; it vanishes in the shadow region. The 
singly- diffracted field is the field diffracted by the tennination of 
the dielectric cover at A (or B). The termination of the ground plane 
at C (or D) yields a doubly-diffracted ray. In this chapter, the 
contribution from each ray will be described in detail. 

A. The Geometrical Optics Field and the Field of the Surface Wave 

Consider a magnetic line source on an infinite, perfectly- 
conducting ground plane covered by a dielectric slab of uniform thick- 
ness as shown in Fig. 2. 



Fig. 2. 
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This configuration simulates the behavior of a narrow slot radiating 
through a dielectric cover. It is know? that the line source excites 
a surface wave which propagates along the dielectric slab. It is a 
straightforward exercise to show that the geometrical optics field Hg 
and the field of the surface wave due to a magnetic line source of 
unit strength are _ 


uEo^r 


-j(kSo- J-) 


Gy. cos %cos + jVe^ - sin^e^sin ?^t 


Hw= - 


^w 

e sin kQ 2 (e^-l) + At(x2e^2 + ^^2)] 


-xx-je|yl 


where ey. = (e'-je")/eo> the relative dielectric constant of the cover 
kg is the free space wave number, k = kg^T^ 

?o * ko^r - sin^0o 


(4) 3 



(5) X 




The value of is determined from the transcendental equation 



In the case where the dielectric slab is lossy, e and k are 
complex, and the root of Eq. (6) is complex. A method^for calculating 
the complex root of a transcendental equation is described in Appendix 
2. It is seen from EqS. (4) and (5) that the important parameters of 
the surface wave 3 and x are complex, so that the surface wave 
attenuates as it propagates. 

Throughout this report a time dependence of e^^^ is assumed and 
suppressed. 
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B. The Fi el d of the Di f f racted Rays 


As noted in the preceding section, the slot excites a surface 
wave which propagates along the dielectric cover (slab). This surface 
wave is diffracted and reflected at the comers A, B shown in Fig. 1, 
so we seek the diffraction and reflection coefficients from the 
canonical problem of the surface wave incident on a dielectric- 
covered perfectly-conducting right-angle bend as shown in Fig. 3a. 
Since an exact solution to this problem is not available, we will 
replace it by an approximate canonical problem, which will be solved 
in two stages. 


In the approximate canonical problem, shown in Fig. 3b, the 
dielectric slab is replaced by the impedance surface described by 
Eq. (1), with k given by Eq. (5). If this impedance boundary is 
positioned at tne top surface of the cover, the field of the surface 
wave external to the cover is unaffected; however, it is apparent that 
the field within the cover is neglected. Thus one would expect the 
approximation to be good only when the ratio K of the power in the 
surface wave field external to the cover to the power within the cover 
is large. It can be shown that 


(7) 


A(k t + % sin^k t) 

A A 

k^ =Jk2 - 02 


It is seen that as t decreases, the ratio increases; and so the approxi- 
mate canonical problem of Fig. 3b is useful only for thin dielectric 
slabs. The behavior of K as a function of the thickness t of the slab 
is shown in Fig. 4, for e^, = 2.56. 


Chu et al [3] have solved the problem shown in Fig. 3b; they give 
expressions for the far-zone diffracted magnetic field and the field 
of the surface wave reflected at the edge. Their diffracted field is 


( 8 ) 


A cos 


Hii =-Hw 


0 I 

xp: 


X + jko cos 0 



-jkoSi 


According to the GTD, 


(9) 


Hil = \ D,, 


-J’Vi 

e 

~^T~ 
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in which is the surface wave field at the edge 
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Fig. 3a. The surface wave incident on a dielectric-covered 
perfectly- con ducting right-angled bend. 



Fig. 3b. Surface impedance model for a dielectric-covered 
perfectly-conducting right-angled bend. 
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SLAB 





Fig. 4. The power ratio of the field of surface wave to 
the field inside the slab. 
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and Dll is the diffraction coefficient. From the two preceding 
equations. 


(10) 


^11 


A cos |- 
X + jkQ cos 0 
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where 

( 11 ) 


A = 


+ X^ 


-SjAp^ + X^ Ij^3-v^|ii/3 X 


^ 2/3 *^0 ® 


j 2 n 


(12) I„ = 


. IT 

2 e sin v0. 


kp sin vir sin 0^ 


. * TT . -1 ^ 

(13) ®s ■ 2" ^ ^ 




The field of the reflected surface wave at the impedance boundary 


IS 


( 14 ) R e 


-j3Sf 


where b is given by Eq. ( 4 ), s Ms the distance along the impedance 
surface measured from the edge, and the reflection coefficient of the 
surface wave is 


( 15 ) R = - 


cos ^65 + it) 
cos ^65 - tt) 
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In treating the radiation from thicker dielectri c covers, where 
the field within the cover cannot be ignored, it is apparent that the 
above approximation is inadequate. In particular, we note that there 
is a significant contribution from the internal field in the form of 
radiation from the vertical end faces of the slab at A-A' and B-B' 
in Fig. 1. 

This radiation from the vertical end faces of the dielectric 
cover is taken into account by using a Kirchhoff-type approximation, 
where the total electric field within the dielectric cover at its end 
face radiates in the presence of a right-angle wedge. Here again the 
dielectric cover is approximated by an impedance surface. The equi- 
valent source of this radiation is the magnetic surface current 


(16) 



_6 

ojee„ 

o r 


(1+R) 


cos kj^(x-t) 
cos k^t 


» 


where k = - 3^ . The magnetic surface current radiates in the 

presence of the wedge structure shown in Fig. 5. R is calculated from 
Eq. (15). To obtain an integral representation for the magnetic field, 
one seeks the Green's function which satisfies the mixed boundary 
condi tions. 



Fig. 5. Line source located at (x = x' e = 0°) on a 

right angled wedge with mixed boundary conditions. 
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9G 

(17a) — -AG=0 , x = 0, y<0 

ax 

and 

aG 

(17b) — = 0 , y = 0, X > 0 

ay : 

This type of the problem can be solved by the introduction of an 
auxiliary function which is a linear combination of the field and its 
cartesian derivatives. The auxiliary function is chosen in such a way 
that it satisfies the wave equation and simple homogeneous boundaries 
conditions. Once the auxiliary function is obtained, the original 
field can be determined by solving a partial differential equation. 

This idea is due to Stoker [4] and Lewy [5] who studied problems in 
water wave theory. The far-zone Green's function, G(s^ ,e;x' ,0) , is 
derived in Appendix 1, where it is shown that 

(18) G(Sj,e;x‘ ,0) 'v- g( 0 ;x',O) 

in which g(e;x',0) is a far-field pattern function which is independent 
of Sj. 

The diffracted field is then 

(19) Hj, = 

where D 12 is the diffraction coefficient associated with the diffracted 
field Appendix 1, 


t cos kv(x'-t) 

(20) 0^2^®) = J3(l+ R) J g(e;x' ,0)dx' 

0 cos kj^t 

This diffracted field contributes a second term to the diffracted 
field previously obtained. Thus the total field diffracted from 
a wedge with a dielectric cover on one surface can be approximated as 


(21) Hj(Sj^,0j^) - H-|^j(sj,ej) + Hj2(si»0i) 
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-jkoSj 

where it is more appropriate to replace 0 by Gj here. With the 
inclusion of the H 12 term, the surface impedance approximation can be 
used to treat the radiation from dielectric covers of increased thick- 
ness, as will be demonstrated in the next section. 

The field diffracted from the termination of the dielectric cover 
is in turn incident on the edge of the ground plane at a distance h 
directly below it. This gives rise to 


-jk^h 


■jk.s 


0^2 


(22) = H^[D^,(|^)+D,2(|I-)]- 








a doubly-diffracted contribution from the edges C and D in Fig. 1. 

Here Dh(0 2 »O) is the hard scalar diffraction coefficient given in 
Reference 6. The factor of % must be introduced at grazing incidence 
( 02 ' = 0) on the edge. This contribution from the lower edges of the 
ground plane has little effect on the pattern in the illuminated 
region, which is dominated by the field directly radiated from the slot 
plus the singly-diffracted contributions from the terminations of the 
dielectric cover. On the other hand, this contribution appears as a 
significant ripple in the pattern of the shadow region below the 
dielectric-covered ground plane. Thus, with reference to Fig. 1, the 
total field for y ^0 . 


(23) 





+ 


+ H2C + 



X > 0 

X < 0 


in which the first subscript 0,1,2 denotes whether the field is inci- 
dent, singly-diffracted or doubly-diffracted and the second subscript 
indicates the edge from which the diffracted ray originates. The ex- 
pression for H when y < 0 is similar. 


A computer program based on Eq. (23) has been developed; it is 
described in Appendix 3. In using this program it should be noted 
that cannot be very close to 1, nor can the thickness of the cover 
t be very close to zero. Thus, the present solution does not reduce 
to the slot in the ground plane without a dielectric cover. The reason 
for this limitation is that the distance from the slot required to 
establish the dominance of the surface wave field increases as 1, 
t -»■ 0. Under these circumstances the field at the termination of the 
cover is not simply that of the surface wave;however, this field can 
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be determined so that the transition to the case of the ground plane 
without a cover can be made. But this must await further work. The 
pattern calculation shown in Fig. 7 is based on a separate solution 
where the surface wave emanating from the slot is replaced by a space 
wave propagating along the surface of the ground plane, and the 
diffraction coefficients for the hard boundary are employed at edges 
A and B. ' . ' • • ' ^ 

III. NUMERICAL RESULTS 

Eq. (23) "is used to calculate the far zone patterns of the don- 
figuratioh shown In-Fig.- 1. In each case the ground plane length L 
is 12.13 inches, the height h is 1.71 inches and its width (In the 
direction perpendicular to the page) is 22 inches. A 0.5 inch by 
0.062 inch slot is positioned 0.1 inch to the right of the center of 
the ground plane, thereby introducing a slight asymmetry in the 
patterns. The axis of the slot is parallel to the edges shown in 
Fig. 1. The precise dimensions are shown in Fig. 6. The slot is 
excited at frequencies so that only the TMq surface wave mode can 
exist. The patterns are calculated in the plane of symmetry normal 
to the axis of the slot for various frequencies and cover thickness . 




0.103" V , 



/‘ / / / / / / / / / / / / // 

713" 

: 1 ,■ 



1 — 1 O 1 ' 

50 " -w 


^ ; ^ 


L- 

-r- = 8.24 AT 8.0 GHz 
Ao , . 

-— = 10.8 AT 10.5 GHz 
.0 


Fig. 6. Experimental dielectric-covered ground plane 
with slot.- • ^ 


The calculated and measured patterns for a slot in the ground 
plane wi -thout a'dielectric cover are Shown in-Fig. 7. As explained at 
the'end of section II, a separate calculation via the geometrical 
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180® F=B.O GHz 


NO DIELECTRIC COVER 


Fig. 7. Pattern of a slot in a perfectly-conducting ground 
plane. Frequency = 8 GHz. 
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MEASURED 

-- CALCULATED 



Fig. 8. Pattern of a slot in a dielectric-covered ground plane. 
Frequency = 8 GHz, t = 0.11 inch, = 2.57. 
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Fig. 9. Pattern of a slot in a dielectric-covered ground 
plane. Frequency = 10.5 GHz. t = 0.11 inch, 

— 2 .56 . 
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MEASURED 

CALCULATED 



Fiq. 10. Pattern of a slot in a dielectric-covered ground 
plane. Frequency = 8 GHz. t =0.1875 inch, 

= 2.56. 
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theory of diffraction is required in this case. Typical calculated 
and measured* patterns for the dielectric-covered slot are shown in 
Figs. 8, 9 and 10; in these figures t^. is the thickest dielectric cover 
in which only the dominant TMq surface wave mode can propagate. In 
the illuminated region, the position and level of the lobes of the two 
patterns are seen to compare in detail, and overall the comparison of 
the two patterns is good. In Fig. 8, where the dielectric cover is 
relatively thin, the end face radiation contribution Hi 2 (oi) of 
Eq. (21) is unimportant and could be neglected. However, for the 
thicker dielectric covers of Figs. 9 and 10, this contribution is 
significant. The importance of H 12 is illustrated in Fig. 11, where 


MEASURED 

CALCULATED 

CALCULATED 



tc=0.4\o 


Fig. 11 


*The measured patterns were provided by Mr. W. F. Croswell of the 
NASA Langley Research Center, Hampton, Va. 
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patterns calculated with and without this contribution are shown in 
the illuminated region. It is seen that this contribution is required 
to bring the calculated pattern into agreement with the measured 
pattern. In this case the slot is positioned in the center of the 
ground plane, so the symmetry of the measured pattern with respect to 
90° is an indication of its accuracy. 

The patterns are very sensitive functions of the parameters, such 
as the thickness of the dielectric cover and the frequency. This is 
shown in Fig. 12, where it is seen that a frequency change of less than 
two percent causes a substantial change in pattern detail. This 
sensitivity to slight changes in the parameters could in part account 
for the differences noted between the calculated and measured patterns. 

The presence of a dielectric cover on the ground plane substan- 
tially increases the pattern level in the shadow region when the 
dielectric is lossless. As the loss tangent, tan 6, of the dielectric 
increases, so does the attenuation of the surface wave incident on the 
termination of the dielectric cover. This in turn weakens the con- 
tributions of the singly- and doubly- diffracted fields. As a result, 
the ripples in the patterns in the illuminated region diminish, and 
the level of radiation in the shadow region is reduced, as shown in 
the calculated patterns of Figs. 13 and 14. In Fig. 14, when tan 6 = 
0.2, the diffracted field components have been so reduced, the pattern 
is essentially that of the di recti y-inci dent (geometrical optics) 
field. 

IV. CONCLUSIONS 

The calculated and measured patterns compare well, which justifies 
the use of the geometrical theory of diffraction in this type of 
problem together with the approximations which have been used to find 
the diffraction coefficient. In particular, the^ method employed to 
extend the approximation of the impedance boundary in treating edge 
diffraction appears to be useful. The patterns of the slot in the 
dielectric covered ground plane are very sensitive functions of the 
parameters, such as the thickness of the dielectric cover and the 
frequency. The presence of the cover substantially increases the 
pattern level in the shadow region when the dielectric is lossless. 

As the loss tangent increases, the ripples in the patterns in the 
illuminated region diminish and the level of radiation in the shadow 
region is reduced. 
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12.0 GHz 



180 ® 

t = 0.190 at 12 GHz 
t = 0.1936 Xjj at 12.2 GHz 


t 


c 


0.4 


Fig. 12. Patterns of a slot in a dielectric-covered ground 
plane: calculated patterns at slightly different 

frequencies. 
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tg=0.4 X 0 


Fig. 13. Patterns of a slot in a dielectric-covered ground plane. 

Frequency = 8 GHz, t = 0.1875 inch, = 2.56 (1-j tan6). 
The values of tan6 are shown on the patterns. 
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Fig. 14. Patterns of a slot in a dielectric-covered ground plane. 

Frequency = 12 GHz, t = 0.1875 inch, = 2.57 (1-j tan6). 
The values of tan6 are shown on the patterns. 
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APPENDIX I 

IKE GREEN'S FUNCTION FOR A RIGHT-ANGLED WEDGE 
WITH AN IMPEDANCE BOUNDARY CONDITION 


Consider a right-angle wedge y = 0, x > 0 and x = 0, y< D, and 
suppose there is a line source located at (p',0) as shown in Fig. 5. 

The Green's function G satisfies the following differential 
equation: 


(A-1) (v2 + 1 ^ 2 ) G(x,y;x',y') = -s(x-x') 6(y-y') 


where 

'x =p cos 0 

* 

y =p sin 0 
and 

|x' = p' cos 0' = p' 

[^y ' = p ' sin 0' = 0 

The boundary conditions of interest here are given by 
dG 

(A-2a) — = 0 ,y=0, x>0 

3y 

9G 

(A-2b) — - XG = 0 , X = 0, y< 0 

ax 


Also, 

(A-2c) G satisfies the Somme rfeld radiation condition, 
and 

(A-2d) G is finite except at (x',y*). 

X is equal to -j kg Zs/Zq, where Zg is the surface impedance given 
by Eq. (1), and in problem of interest here, is such that RpX > 0. 

Our treatment of this problem closely parallels that of Karp [7], 
except for the manner in which the unknown constant A, which will be 
encountered later, is found. 
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Let us make a transformation 


(A-3) V(x,y) = (|- - a) G(x,y) ; 

9X 

then Y must satisfy the following differential equation: 
(A-4) (v2 + Rq2) v(x,y) = (~+ a) 6(x-x') 6(y-y') 

Using the relationship ( 9 / 3 x) 6 (x-x') = -( 3 / 9 x') 6 (x-x') , 


(A-5) (v2 + ko^) V(x,y) = (-^ + a) 6(x-x') &{y-y') 

dX 

The boundary conditions for V(x,y) are then 
9V 

(A-6a) — = 0 , y = 0, X > 0 

ay 


(A-6b) V = 0 , X = 0. y< 0 


The function T(x,y) need not be finite at the origin since it 
involves a differentiation of G. The solution to Eq. (A-5) can be 
obtained by first solving for the Green's function GQ(x,y;x' ,y' ), which 
satisfies the differential equation 


(A-7) (v2 + ko2) Go(x,y;x' ,y') = - 6(x-x') siy-y') 


together with the same boundary conditions 
aG 

(A-8a) = 0 >y = 0, x>0 

ay 

(A-8b) Gq = 0 , X = 0, y< 0 


and the radiation condition. 


24 



It can be shown that 


(A-9) Go(x,y;x* ,y') = ^ cos 0 cos e' Jp ^ (kp^ ) - 

n=0 


2n+l'"P< 


H2n+i(kp>). 


( 2 ) 

where i s Bessel function of order' v and H - * is ;a Hankel functi on of 
second kind of order v. v 

By applying the operator -(9/^x'' + x) to GQ(xjy;‘x‘ ,y ' ) , we- can ' 
obtain the following solution for V(x,y): 


(A-10) V(x,y;x',y') = + A) GQ(x,y;x' ,y ’ ) 

The complete solution for iT{x,y) is 


(A-11) ' V(x,y;x'-,y') = %(x^ 


-'■...1 ••• ! : .+, -I an -cos e' (kp ) ;• • ' - '* 

n=0 ^ 

The last summation contains the-cl ass of admssible singular solutions ■ 
'which satisfy the homogeneous equation for V. Since V cannot be too 
singular at the edge, the only admissible function of the singular class 
is the first term; hence, 


(A-12) 


dn C 1 


= 0 , n > 1 


A particular solution of (A-3) is 


(A-13) Gp(x,y) = -e^"" J e"^^ V{§,y) d§ . 

X 


The existence of the integral is guaranteed since Re X > 0. 
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The function Gp(x,y) does not satisfy all the conditions of the 
. problem, for it is hot a wave function on the negative x-axis, since its 
y derivative is not continuous there. This defect can be remedied by 
adding a complementary solution for 6, obtained from (A-3), to (A-13); 
thus we obtain 

/ Xx+ j y 

026 , y< 0 

(A-14) G(x,y;x',yV) = G„ + ^ 

1 0 . y > 0 

-Cj and C 2 are then determined by imposing the boundary conditions 
(A-15a) G(x,0'^) - G(x,0") = 0 , x< 0 


and 

aG(x,0+) aG(x,0") 

(A-15b) 0 , x< 0 . 

ay ay 

Numerically it is very difficult to obtain Cj and C 2 . This difficulty 
can be avoided if we choose a suitable function [7], H(x,y), which has 
the same singular property as H ‘2) (k^p) and is a liomogeneous solution 
for V. 1/30 

Let us consider the Green's function N which satisfies 


(A-16) (v2 + ko^) N(x,y;x',y') = - 6(x-x') 6(y-y') 


with the boundary conditions 

aN 

(A-17a) — =0 y=0, x>0 

3y 


aN 

(A-17b) — = 0 X = 0, y< 0 

ax 

together with the 


(A-17c) Radiation Condition. 


26 



It can be shown that 


(A-18) 


N(x,y;x',y') = - ^ I "T® 

3n=0 " ■ 


cos 




j^(kp,) H|;^(kp;) 
3 3 


where is the Neumann's constant 



n = 0 

n = 1,2, ••• . 


; 


(A-19) 

(A-20) (v2 + ko^) H = 0 . 

H(x,y) is thus a solution of the homogeneous wave equation with the 


boundary conditions 


(A-21a) = 0 

3y 

II 

o 

X 

V 

o 

(A-21b) H = 0 

X = 0, y< f) 

together with 



(A-21c) the Radiation Condition. 

Note that the leading term of 8N/3x has the same singular property as 
The complete solution for ¥ is then 
(A-22) V(x,y;x',y') = -(^+ A) Gp(x,y;x' ,y ') + A H(x,y;x',y') 

where A is a constant to be determined. 
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That G satisfies the wave equation can be shown from the fact that 
Gx = V + X G 
G^x = V + A V + G 

also Vyy = -ITjj - X) 6(5-x') 6(y-y') 

®xx * Sy ' '■ 

_ _ 

V + X V + x2 G - e J e V (§,y)d§ + k 2 G . 

X o 

Integrating by parts twice in the above integral, 

3 A(x-x-) r 

®xx ^yy "^o^ ^ ® J «5(§-x') d§ 

. X(x-x') 

= -(^ + ^) e u(x-x') 6(y-y') 

X(x-x') „ 

= -e ^ u(x-x') 6(y-y') 


= 6(x-x') 6(y-y') 


where u(x-x') is a unit step function. 

In the far zone, we observe that Gq, N, G and ¥ have the asymptotic 

form 

(A-23a) Gq - — - — go(0;x',y') 

y/p 

^-jkoP 

(A-23b) N ^2 n(0;x',y') 

(A-23c) G - ^ g(9;x',y') 

y/p 

and 
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(A-23d) V u(0;x',y') . 

4 " 


Since the far-zone field is of interest, only those terms which are 
ofO(l//p0 are retained, the operator 3/9x can be approximated by 
-j k cos" 0; thus we obtain 


(A-24) 


G(x,y;x',y■)<i^ 


[(A-1)-^ - x] Gq + j k(j cos e A N 
j kQ cos 0 - X 


Hwang [*] determines the constant A by transforming the function Gq 
and N into the integral representations which are then evaluated asymp- 
totically by the modified Pauli-Clemmow method of steepest descent. From 
the asymptotic solution, we can identify the ray-optical behavior of the 
field. The pole contributions give rise to the geometrical optics field 
which is the sum of the incident field and the reflected field. From 
this we obtain 


(A-25) A = 1 . 

Thus, 

j k- cos 0 N - X Gq 

(A-26) G(x,y;x',y') ^ — — 

j kp cos 0 - X 

in the far zone, where N and Gq are given by Eqs. (23a) and (23b). 


*Hwang, Y.M., "Electromagnetic and Scalar Diffraction by a Right-angled 
Wedge with a Uniform Surface Impedance", to appear. 
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APPENDIX II 

ON CALCULATING THE COMPLEX ROOT OF A TRANSCENDENTAL EQUATION 


The problem of finding the complex roots of a transcendental 
equation is usually reduced to the problem of finding the roots of two 
simultaneous nonlinear equations. The two nonlinear equations occur 
because the real and imaginary part of a complex function must also 
equal zero when the complex function f(z) = 0. The method of finding 
the roots of simultaneous, nonlinear equations has been discussed in 
many books; but sometimes it is tedious to reduce a complex function to 
two real nonlinear functions. In this appendix we point out that 
Newton's method for the real roots of a real function can be extended 
to find the complex roots of a complex function. 

Suppose that we seek a solution Zq of a function f(Z,e) such that 


(A-27) f(Zo,e) = 0 , 

where e is a comple^ parameter. 

Using Taylor's expansion, 

(A-28) f(Zo.e) = f(Zol.e) + ^'(^ 01 .^) 


where Z^j is an initial estimate of the root, which is assumed to be 
close to Z . For example, Zq^ could be found by solving a related 
problem exSctly. 

Let 


(A-29) aZj Z.q - Zq^ 

From Eqs. (A-27) and (A-28)', 

f(7-oi,^) 

(A-30) aZ. = - ^ - aZ 


Next, replace Z by Z^^ = Z^^ + aZ^^, and repeat the procedure; in 
this way we can generate successive approximation, and after N iterations. 
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(A-31) 


Z 

0 



N 

I 

i=l 


az 


1 a 


» 


where 

(A-32) AZia 


f (Zoi »e ) 


When Zqj^ is chosen sufficiently close to the root Zq and 
i 0, the process will converge, i.e., < I'AZ^gl. The 

procedure can then be terminated when |AZ]aj^< S', where 5 is a 
small number. 


f'(Zoi,e) 
i terati ve 
sped fied 
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APPENDIX III 

THE FORTRAN IV PROGRAM FOR CALCULATING 
THE FAR-ZONE RADIATION PATTERN- 


C PATTERN CALCHLATUl.xi '.'E S LU "i KAO I^TOR IN a“T)I 4L EC’Ir I C-C"A T^Ti U l'~p1' ATE 

_C THIS PRCGR AI- IS r , nlli: u-gY Fl'K 0 < n < 0 L Hi o 

C LA>OAO = FKFF " SP/''r.E S’ A VF ' LEN (-f H S ■ 

r. THIS PRDGRAFi Cl'-''inT AF liSFO '.'.'HtAi F S< = 1 . U AHn/nQ l) = 0.0 

C 

_C FSK - RELATIVE niFLFCT PIC Cll hSTAMT ^ 

C FSO - FREE SPATE f) I F LF CTR I C Cll'-i ST AW f 

_C FL - niFLFr. TRi r. C O' -ST AHT OF UT S sV S L A H 

C n - THICKVESS nF“,ilETECTRIC O'aHmt; Tn IA'Chfs 

_C F - FRFOU-I'ICY IT' CPS 

X B - FRFF-SPACF '.;AVF i-i.k.kFR 

_C BK - DIE L ECTRIC '-■iav E ^'liHAF R Ff iR LO S SLE SS SLAJ^ 

C VI - THICKmFSS cF TRlirCATFl) GRfHii'iD PLANE Tm I.'X'hFSIHoT iMCLI'l'iiJG 

_C the thick ' E S .S_! 1F_ S_UiB J ^ 

C HI - TIIOTH of niELFCTP.lc COATI'h; FORM THf" SniiRCF TO THE RIGHT 

_C CnR^■ER IM INCHES 

C XL - TOTAL HIPHT OF DIELECTRIC COATING I iM IMCHES 

_C XLT - LOSS TAiVG-MT (iE LC'SS Y SLAB ^ 

C T - COHVERSir,, EACTL'Hll fl'iCH = 0 .025 AME T FR S ) 

_C Cll - SlIRFA CE-!'' AVE ATTEi^'IIAT lO 'V C l.iHS I AiX[ 

C CBE.TA- SORFACE-mavE PROP i->, AT foFi CONSTANT 

C ■ HGO. - DIRECT RAOIATIQM rlE L D 

C FJIU ■ OIFFRACTiniv FIELt'MKIGHT SIDE) 

C HD2 - OIFER A CTin i -! F I E L O ( LEFT SIDE ) 

C HY - slab R ADI ATKlf',' FIELD 

C ■ HT - total RADIATIOh F I E LD 

C FN - NORMALIZED CUNSTAMT CJE THE TOTAL RAOIATIDiv FIELD IN DK AT 

_C 90-degree S 

C HTT - TOTAL RAD I AT If!'", FIELD INCLUDING SLAE RADlATIuN 

_C DRHTT-TOTAL P'AITIATION FI S lD IN DB( W1 TH SLAB RADI ATI UN) 

C PHD - ASPECT angle In DE'gREES ’ ' ’ 

_c : 

COMMON BJ( 150 ) ,RY ( 150) 

n IMENS ION A1IH( IRG ) ,AUHS ( IBQ) ■ ■ 

COMPLEX HGO ,TH£TS,B1 , 8 2 ,B3 , A 1 , A 2, A , HD,H T, CCS'-' *HR,HLf HDl ,H02,HY 

complex hot 

complex CN f CO ,C2 .HTT 

CiDMPLEX HYF ,HYR ,GF ,GB 

COMPLFX CCSV.'l ,CCS'H2 

COMPLEX HDID ,H0 7n ,HR1 ,HR11,HL1,HL11 

COMPLEX HP, HO 

COHPLE X OHl ,0H2 ,0H 3 , DH4 , i)M5.,H 

COMPLEX HO in ,HD1 1 ,HD20 ,HD21 ,DD ,DVDR ,DV DL 

COMPLEX CSH,CI 1 .CI2.X 

complex EL,Z ,ZK,C0,CBETA,CES1L,CESR,TS,C1 

Pl=3. 141592 , ■ 

1212 R6An(5,llll .EM0 = 1112)ESR,F ,FN,0,XL,W1,V1,XLI'’ 

nil F OR MATI F4 .2 , E1D.3 , F6.2,AF7.4,E9.? ) 

U2 = XL-W1 

URITE (6 , 1313 )F , 0 ,H 1 .iV 2 

1313 FORHAT ( IHl , ' FREOUEi'!CY= ' ,E 10.3 , >HZ SLAB TH I CX ME S S= ' , F F. A ,'I;'’CH 


1 Wl=< ,EB.A ,5X, ? = < ,EB.4 1 


€ 

WRITE (fc,BfTHB )XLT- 

ARAB FORMAT!' LOSS T AM GENT = ' , E 15 .A / ) 


WR ITE ( 6, 1444 ) 

1444 FORMAT!' PHD - ASPECT ANGLE IN DEGREES'/' ITBH U - TOTAL RADIATIf! 

4 

IN FIELD IM DR'//) 
B = 2.0=:=PI*F/3 .OEOB 

3 

FSU = A.A5E-12 
T = 2.54E-2 

2 

i.il=Wlj;=T 

W2=H2X--T 
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I)=|)*T 

yjj:Vl±T 

V = (FSR-1 .0)<=( r 

CALL FhP(F ,F<;R,n ,y ,h S 1 L . -F T A ,U , 1 ) 

Xi)=FSlL-n 

aa=fs^-.=f so : 

P=XLT«AA 

LLf CALEL.X (.AA.ltP ) 

call HSl’ oTn,F,£LtXir,7. ,ifcp) 

7K=( (?.=:=£ p.^F )g*?)*^.~Pl»l.F-7»£L 

CH£TA = CSORT (Z<-( 7.’rZ / U>=>:'< ) ) ) 

Cn = CSOR T ( ZK- B*B - ( Z *Z/( U »C ) ) ) ^ 

cesi L =Z/I1 

CFSR=F L/E Sf) ^ 

CHBTA = Cni\'JC;(CRETA) 

CU = CU'^i JG( Cl) ) 

CESI L = COi'!JG (CFSI L ) 

■ CESR=CO'':J Gir.FSR ) 

CALL HTSO'j ,Ol , B, TS, lER ) 

C C? IS THE BEFLFCTIOf'' Cf lE FFlCI Ej^T 

C2--CCnS( ( TS + Pl ) /3. )/CCOs’( (TS-PI ) /3 .) 

1^0 : 

PHO i-l.O 

1 PHD = PHO 2. 0 

PH = PHr)=?n.ni7AS 

1 = 1+1 . 

X = B*CSORT ( CESR-CnS( PH ) s=Ci!S ( PH ) ) *0 

HGU = -2.=:=Pp.-F*CESR’:^FSti*Sr;(PHI^CFXP(CAPLX( .O.-PI/A'. I l/( (CFSR»SIN(PH 

1 )-'=CCUS( X )-CrPLX( .0, 1 . ) =;‘CSC'kT (CESR-COS( PH )5;=cn S( PH ) )*CSI H( X ) )=:= 

IS ORT ( 2 .»P 1 »R ) ) 

CSW=-2.*PI -F=:=CESR*CESR-ESO*CU=i=CU*CE SIL/( C8 E 1 A*CS I M( C ESI L*0)*(CESK» 

1 ( Z K- 3 ) + CU -lA » ( r.H =i CII »CE SR »CE SR ±CE$1L=P CESI LI ) ) 

CCSW 1 = CSH=.K6XP ( C'-'PLX ( . 0, 1 . )«CBET A*W 1 ) 

CCSH2 = CSH»CE XP ( C-IPLX( .0 , I. )xCBETA*h 2) 

Cl 1 = 2.*CEXP (0-1 PLX t .0,-PI /6. ) )*CSli'! < TS/3. )/( B*SIN( PI / 3. I-CSIi'M TS ) ) 

C 1 2 = 2. -CFXPI CM PLX ( .0,-P 1/3. ) )^CSIM(2.X‘TS/ :-.)/( s I 'V ( 2. -•■■=P 1/ 3. ) =:= 

ICSIM (TS) ) 

Cl = CSORT( B--:--R+CU*CU ) 

A=CMPLX ( .0,-4 . )*Cl/( a-'PLX( .0,3 .) -C 1 *CI 1- SOR T (?.)*( Cl l*Cll-CI2*R- 

1CEXP(CMPLX( .0,2.4-p t/3 . ) ) ) ) 

HO = A -SORT (2 ./ t PI -R ) ) -CFXF { CH P LX ( . 0, -5 . *P I / 1 2 . ) ) /( -Cu-CMPL X( .0 , b) - 

IS IN (PH)) ' ^ 

El=R=;='.-ii 

F2 = B»!--'2 ; : 

IF (PH0.LT.1P0. )G0 TO 4441 

0H= ( PHO-3 AO. )».0l74'j ; 

G!l Tl) 4442 

444 1 f)H= PH 

444 2 HD 1 =CCSWl'=HD«CnS( ( P I / 2 .O+OH ) /3 . 0 ) *C6 X P ( CMP L X ( 0. 0, -E 1 *C OS ( PH ) )T' 

H0 2 =CCSl.'2»Hn--;=CnS ( ( 3 .0«P l/2. O-PH ) /3 . 0)4CEXP ( C mPL X( O.O, F2-~:-CnS( PM) ) i. 

XM=1,5 

! XLGDA = 3 .OFOP/F 

nn=A=;=SORT(2 ./( PI -p ) l^CEXPICi-iPLXI .0,-5.*P I / !?. ) ) / (-Cu + CHPL X( .0,8) ) 

« HO io=ccs^-' 1 4-no ^ 

Hi)20=rxsi'2*(in 

! IF (PHn.LT.90. )Gn TO 661 

pHV=4S0.-PHn 

t Gf) T(l 662 

661 PHV=90.-PH0 

j 662 ph;-;=-oo. + pho 

Yvn= (D+Vl ) /XLHDA 

£ CALL OFRCF (nVDR , X/M , Y VO , PHV , X L W A ) 
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CALL OFRCF n'^riL.X'', YvT- , ■" ■ - ■ -- - — 

n'./L)R = r.iiA!j(;{|iV'!)i^ ) 

in/o L = 0 )'v jc (’n'vh. L I " 

Hi)Il.=Mino::=l lvi^H=^r.P XP [ C'-PLX ( ( ) -t-V 1 )i= ( l.O-t-Sb'K PH ) ) -" CnS ( F H ) ) 

1 ) ) /suPT ('1+y 1 ) 

H fi? l = Hi) ?0»:=iiynL--:=r. F XP( Q - PLX ( o .o, B* ( ( ) +V 1 jfj UO+SLr'( >'H ) ) +!’2--:=CnS( PH) ) 

n ) /Fi.'RT ( vi+tn 

r, THIS Sl) PK ^'niTI- iF COMPII T FS THE SL AP RADIATl OfM FIELD ^ 

PHF=PHn+PD.O 

\-: = 2 .O-'-l I 

CALL HSLAP (O ,v 1 ,F SR ,CU ,CHE T A , CF SI L . CC S"' 1 ,PHF ,hyF ) 

- hyf =hy f =!=cf xp ( c-i PL x ( n.D, i«c u s( ph > > > 

PHH=270. O-PHO 

V.' = 2 .0».'2 

CALL hslar (n ,v I ,f ,!•' ,CFSR ,CU ,'cbET A,CE sl L . CCSW 2 , PHH ,ri Y y ) 

HYK = HYR=:^CE xp ( C‘'| PLX ( 0 .(),B»'2»CU S( PH ) ) ) 

hy=hyf+hyb 

HY = HY»( 1 . 0 -C2 ) 

IF (PHI1.lt. 90.0 )G;i TO 2221 

IF ( PHO . LT . 1 p. 0 ._n )_0f I 2 2_22 

IF (PHn.LT.r7n.0)0ii TO 2223 
HT = HU I +H0 1 1+HD2 1 
•00 TU 222A 

2221- HT = HG0 + Hlll+HO2+Mni 1 
c,n TU 2224 

2222 HT = HGfi+HO l•^Hn2+'-l^2 1 
GO TO 2224 

2223 HT=HQ2-i-Hr)21-fHnir 

2224 HTT=HT+HY 

AHT = CABS(HT ) ; 

AHTT = CARS(HTT) 

nBHT = 20 . yalOGIO ( AHT 1 

riBHTT = 20 ,0=!'-ALnG10( AHTT ) 

AUH( 1 I =OHHT 

AUHS( I )=()BHTT 

IF (PHP . L T. 35 9.0 )00 TO 1 

FNl=(AUHSl4!3)-t-AliHS(46) )/2. 

Y2=FiM-FM 1 

ANGLE=-1. 

003333 I = 1,1 BO ^ 

AMGLc =ANC,lE +2 .0 

f1BHT = AUH ( 1 ) +Y2 

npHTT=AUHS( I ) +Y2 

WRITE (A, 4 ^ 4 ) AM GLF ,nRHT>OBH TT 

444 FORMAT!' p’hO = ' ,F 9 .2 , 2X , T'EHT = • ,F 10 . 3. 2X , ' ORhTi = ' , Fl'o . J) 

; 3333 COinTIMUF ^ 

GO TU 1212 ■ ' 

1 112 STUP 

EMI) 

SOPR OOT lis'F HT S (0 ,CU,R t T S, IE R ) 

COMPLEX CO ,fs 7f.SC i TSI fbcu ,TSR 

7 CY=o/n ^ 

CALL TH£T(CY,CI) 

6 TS0 = CHPLX(R.141SQ?/2. . -C 1 1 ^ 

n£L=l.E-3 

5 A = ATO A G( r.i I )/ R F Al 1 CO ) 

A A = A B S { A ) 

3 TSI^^TSU 

IF ( AA.LE .OE L )G0 TO 10 

J L lK;LL=CM^±)a.P.t..^LiXiil.RO:LlCL'.J*i%ll£L/J1i^^ 

CALL CRTSIOCO , K,T SI ,TSH . ItR ) 

; IF ( IFR .EO. 1 )G0 TO 3 
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TSl =TS« 

neL= QFL -H .F-‘S 

IF '( A/i.GT .OE L )('^n TU 1 

10 CALL CKTS<au>^,TSl ,TS.1FR) 

3 RETUKiV 

FJ^) , 

"sDHRunT iMg CFTS( rjj’, =.,t’?'^tsVier ) '" 

Cnf-iP LFX F , CU «T S .nT?l! ,7S ! S iJ SJ ^ r 

F ( CM ,R,TC) = (CSORT ( 1 .+Cl’«CU/( F*rl) )-CsIn ( TC ) ) / CCnS MC ) 


M = 0 


DTSO = F ( CU ,R ,TSn ) 
Tl = CARS(nTSll ) 



TSI=T.SO+OTSr' 



1 _ 

A=CABS( CSORT ( 1 . +Cl!«ai /( ) )-CSIN( TS I ) ) 


PTSI =F ( Cn ,R ,TS I ) 
T? = r.AB.s(nTSi ) 

IF ( A .LT . 1 .0-5 ) (i!) TO 2 
IF (T2.GT.Tl)r.ri m 3 

T1 = T2 

TS1=TSI+I)TSI 

IF (•■'!. GE .30 )GO TO 5 



GO TU 1 



2 

I£R=0 


TS=TS1 
GO TU 7 


3 

1ER = 1 

WRITE (6,4 ) 



4 

FDRiiATI ■ THIS ITERATION OOES NOT CONVERGE') 
GO TU 7 



5 

IER = 1 

WRITE (6,6) 



6 

FORMATC THIS COM. PU T AT lOW DO E S NOT CONVERGE IN JO ITERATIONS') 



7 

RETURN 


EM) 

SURRUllTIiME HS'-) Pin ,F ,E L, XI) ,Z , lER ) 

COMPLEX EL,Z ,RZ,ZI ,ZO,EP 
Zn=CMPLXIXO, .0) 

nEL=l.E-5 
AIE L=A TMAG (E L ) 

REL = RFAL(EL) 

A = ARS ( A IE L )/RE L 



Z 1 =zu 

IF ( A.LE .OEDGO TO 10 

\ 


1 

EP = C0PLX(RE L',-RFL<=nE L) 

CALL CRSH Pin ,F,EP,Z[,RZ,IER) 


IF ( lER.EO.l )GU TO 3 
ZI=RZ 

0EL=nFL+l.E-5 

IF ( A .GT .OFUGO TO 1 


10 

CALL CRSWP(I),F,EL,ZI ,Z ,IER) 


7 

3 

RETD RN 




Ei'G; 


e 


SORROOTIMF CRSWP(I) ,H,EL ,Z0 ,RZ ,1ER) 




COi.PLFX Zn,HZ ,EL,Z,E,G,FP,GP,DZO,nFG,Zl , DZ 1 ,ZK,E1,ZZ 


5 


F (Z)=Z*CSIN(Z ) 




G(E L ,FSrj ,ZK ,R ,0 ,Z ) =F L« CSORT ( (ZK-R«h )*n=?l)-Z#Z )-r.C(IS(Z ) / ESO 


i 


FP( Z ) =Z-CCOS ( Z )+CS (Ml Z ) 




GP (F L,t SO ,ZK ,R ,n ,Z) =-( FL/E SO )* (Z-CQISI Z ) + ( ( Z K - R4 (s ) =; n*n-Z *Z )* 


i 


ICS IN (Z ) ) / CSnPT ( ( ZK-R*H) *l)*r)-Z«Z ) 




ESO=H.85F-l? 


t 


PI =3. 1^15 0 
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h = 2.*PI#H/3 08 

ZK= ( (2.«PI =i=H )«=?P )=i=4.=sPI *1 .E-7*EL 

r.i = 0 . 

0 70=- (E ( 70 )-r,{EL,ESQ ,7K ,R,t) ,70 )» /{FP ( 70 ) -GP ( F L , ESP , 7K, B, n, 70) ) . 

Tl = CA8S(l)7n ) 
7I=7U+D7n 


1 

DFG = F (71 )-G(E L,FS0,7K,3,L),71 ) 
A=CABS(OFG) 

D7I=-(F ( 71 )-G( FL,ESn,7K,B,D,7I ) )/(FP(7I )-GP( FL,FS0,7 K,R,0,7I ) ) 
T2 = CARS(P7I ) 

- 


IF (A.LT.l.E-5 )GP TO 2 
IF (T2.GT.T1 )G0 TP 3 

Tl=T-2 

71=71+071 

N =N + 1 

IF (N.GE .30)GP TP 5 


2 

GO TO 1 

ier'=o 

R7=ZI 
GO TO 7 


3 

IER = 1 

WRITE (6,4 ) • 


A 

■FORMAT(' THIS ITFRATION DOES NOT CONVERGE IN SEARCHING FOR COMPLEX 
1 ROOT OF SURFACE HAVE POLE') 


5 

GO TO 7 
IER = 1 


6 

WRITE (6,6) 

FORMAT)' THIS COMPUTATION DOES NOT CONVERGE IN :-0 ITERATIONS IN SE 


7 

lARCHING FOR COMPLEX ROOT (iF SURFACE WAVE POLE') 
RETURN 

6N0 

SUBROUTINE SWP (7 ,ESR,0,V ,£SIL,HETA,U ,IER) 


C 

C 

THIS ROUTINE IS TO FIND THE ROOT OF THE POLE OF SURFACE WAVE BY 
USING THE BISECTION METHOC. EE SURE 10 CHECK THE FIRST Tl-P 


C 

APPROXIMATIONS XI AND X2 SUCH THAI F ( X 1 ) 4F ( X2 ) IS NEC^ATIVE 
FIX) = TAM(X) - ESR *SORT ( (V/IX*X ))-1.0) 

R = 2.0<=3.1A159«7/3.0E08 
RK=SORT (E SR ) =5=8 

IF (O.LT.2.5AE-5 )GU TO 1000 
N = 0 

X1 = S0RT(V)>!‘1.E-A 
TV=1.5708 

SV=SORT (V) 

IF ( SV .GF .TV )G0 TO 1111 

X2 = SrjRT (V )<=.9999 
GO TO 992 


1111 

992 

X2=TV-0.0001 
Y1=F (XI) 


^ 1 

Y2 = F(X2) 

X = (X1+X2 )/?.0 

7 


N = N + 1 

IF (N.GT.30) GO TO 2 

C 

100 

IF (F ( XI) «F (X ) ) 100,10 1, )02 
X2 = X 

5 


lF((X2-Xl).LT.1.0E-5) (0 TO 101 ^ 

GO TO 1 

4 

102 

XI .= X 

IF((X2-Xl).LT.l.nE-5) G(.i TO 101 

3 

101 

GO TO 1 
ESI L = X/I) 

2 


BETA = SORT IBK=i==x2-ESl L*-2 ) 
U = SORT{RETA*?;s2-R-*2 ) 
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r-u Tu 1001 

'1000 ll= (FS‘^-l 0)*H* flg()/ESK ^ 

wETA = Sl'«T Ul*'l + t<*r, ) 

FSII. = f>»S0m(F.SR-1 .01 

1001 1EK=0 

Oil TU 3 ^ ■ 

2 Via = 1 

PRINT S ^ 

A FnRHATl • ','FAILFf) TO CiiMVEROE IM 30 frERATIONS IN SEARCHING FOR. 
ITHE ROOT OF THE POLE OF SURFACE W AVi » 1 

3 RETURN 

ENO 

.Subroutine tnetip,o 

C THIS ROUTINE IS TO COMPUTE A COM S T AM T BY F I NQl NG A ROOT OF A 

C POLYNOMIAL BY USING THF SECAND METHOD IN ORDER TO COMPUTE THE 

C CHU-KfUIYOIIMJ I AN <S CUE FF IClEiMT 

FIX) = SINHIX l-P 

X = 0.0 

XI = P 

Y = FIX) 

Y1 = F IXl) 

nn 1 I ■ = 2,20 

X2 = IX«FIX1)-X1=!=F|X))/(F|X1)-FIX)) 

Y g F 1X2 ) 

IFIAHSIX2-XD.LT.1.E-5) GO TO 2 
X = XI 

1 XI = X2 

MRTTF I 6. IS) . 

6 F()RMAT(/36H0FAILED TO CONVERGE IN 20 ITERATIONS) 

GO TO 3 

2 C = X2 

3 RETURN 

end 

. SUBROUTINE DFRCF (D ,XN ,Y ,P ,XLMD A) 

C THIS ROUTINE IS TO COMPUTE THE DIFFRACTION COEFFICIENT 

complex F lJ,F2J,0,R,Tl,ni ,T2.D2 ,0,F1 ,F 2.S ' 

PI=3. IA15927 

BR = B»0-.017A5329 

ARG1=IPI+RR )/I2.0*XN ) 

ARG2=IPI-BR)/I2.0»XN) 

CX1=C0SI ARGl ) 

CX2 = CnSIARG2 ) ^ 

SX1 = SIN| ARGl ) .. 

SX2 = S1N I ARG2 ) 

X1 = IBR+P1 )/I2.0’S‘XN<=PI ) 

N1 = X1 

g 1=X1-N1 ^ ^ ^ 

IF IEl.GT.0.5 ) iMl=Nl + l 


IF (E l.LT .-0.5 ) N1=N1-1 
EM1=EL0AT INI) 


7 

X2=IBR-PI )/I?.0«XN>t=PI ) 
N2=X2 . . 


< 

E2=X2-N2 

IF IE2.GT.0.5 ) N2=N2+I 


S 

1FIE2.lt .-0.5 ) N2=N2-1 
EN2=ELDAT IN2 ) 


4 

A1=1.0+CDSI-«R + 2.0=!‘XN«P1vPM1) 

A2=l .O+COSI-BR+2 .0«XN*PI =f=FN2) 


i 

SA1 = S0RT 1 A1 ) 
SA2=SORTIA2) 


2 

XX=I SORT 12 .OA-Pl #Y ) )«SA 1 
YY=ISORTI 2.0*P1 «Y ) )«SA2 
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xxs=xx-xx 

YYS = YYS=YY 

pn = S(.'RT(PI/2.0) ^ ^ 

. CALL r.SCCl ,S1. ,XXS ) 

cci=o.ij-a 

.. . .SS1 = 0.5-S1 ^ ^ 

F1J= PO=^OiPLX( SSI ,CC1 ) 

CALL CS(C?,S?,YVS ) 

CC2=0.5-C? 

SS2=0.5-S2 

F2J= Pu=i=CMPLX( SS< ,CC2) 

P--P1 /<♦ ,n 

a=C«PLX (0.0 ,P ) 

R = CFXP( (.) ) 

S=(-0.25 yksSOKT ( XL'-'OA ) )/( PI s=X^ ) 

Tl = S=!=? .n»CEXP( C<PLX ( 0.0 ,XX5 ) ) 

IU=Tl=i=F 1J=:=CXIY( xx/SXl ) 

^ T2 = SX=?.QX=ceXP( OYPLX(0.0 ,YYS)> 

0 2=T2X=F2 J-CX2*( YY/SX2 ) 

n= ni-n)2 

RETURN 

END' 

SllBRUliTINE HS LAP (0 ,V 1 ,F ,1- ,CESR .CU ,CBE TA , C ES I L t CCS" » PriO,H 1 ) 

C fHIS ROUTINE is TO COMPUTE THE SLAP RADIATION FIELD 

COMPLEX CCSU.HYl ,HT,CESR,CU,CHETA,CESiL.X‘kX ,XKXC),z< 

PI = 3 .lA 159 ; 

ESU = H.R5E-12 

XLHDA=3 .OFQP/F ^ 

R=2.0*P1/XLMDA 

ZK = CE SR«H»q ^ 

XKX=CSCIRT ( ZK-CBETA0C8E TA ) 

. XKxn=XKx*n ^ 

,HT=CMPLX(0. 0,0.0) 

. ' H=n/ (0.05»xLMnA ) ^ 

on 222 I = 1,IM 

M = 2»l-1 

ni=o=N'./(2.*iM ) 

. CALL CASER(D ,F ,0 ,V 1 tCESR ,CU ,CBET A,Ce SI L ,X L NOA , PI , PH |;,HY1 ) 

HT=HT+HYl*CCnS( XKXYID-DI ) ) 

222 COnTIMUE 

HT=CUNJG(HT) 

HT=HT*0=i=CMPLX( .0,-1. )*CHET A»CCSW /( 1M*CC US ( XKXD) ) 

RETURN 

END 

SUBRUOT INE CASER (O ,F ,H ,V1 fCESR ,CU ,CME TAtCESI L.XLMOA, DI ,PHn,H ) 

COMPLEX riHl,nHYTnH2,H,HY,hni,HI)2,FG,OB,CESR,CO,CHFTA,CFS IL.UK 

complex C 

PI=3.IA159 ^ ^ 

BK=2.0*P1/XLMDA 

IF (PHD ,GT. 270. )PHn = PHI)-360. 

PH=PHD<=.017A5 

2 UK = CLJ/BK 

ARG .= BK.<=UI • 

« CALL GREENF (PHD , ARG ,UK , X LMD A.FG) 

HY=FG 

5 C = CEXP( CMPLX( .0,-Pl /A. ) )/SORT ( 2 .»PI»PK) ; 1 

XN=1.5 

1 Y=(l)-l-Vl)-ni ^ ■ - 

YO=Y/XLMi)A 

2 PHl = lR0.n-PHn : ■ ' ■ ' ' 

CALL DFRCF (DHl ,XY ,Y0 ,PH1 ,XLNDA) 

» HI)l=C»llHT»CEXP(O'-PLX(0.0,(HK»(D-*-Vl) »C nS(PH )-H.K=XY )1)/ S"RT(Y) 
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CALL nPRCF (HHY ,X'' ,vn ,O.C,XLf'PA) 

Y? = (. 1 .) +y_l )_/ X Ljt\Ll.A 

CALL "ni^RCF (riH2,x - ,Y2 ,Ph'i ,XLf-:i>A) 

Hn2 = C*QHY=:=|lH?x=CF XP (C(. PLX ( .0, ( -BK<=‘(I)4-V 1 )-bKYY ) ) ) /SOR 1 1 Y * I l>i- 1 ) ) 

IF ( PHI' ,LT .0 . )G(1 T(! lli:- 

1F(PHH .LT, l_fl 0^ ) PM TO J. _1J_1 ^ 

H = HY+H|)2 

r.ii TO 2 

1U3 H=Hni 

CO TO 2 

nil H = HY + H!) 1+HD2 

2 CONTI MU F , 

RFTORW 

FNO 

SUBRUUTIMF CRFEA'F (PHD,ARr,,UK,XLHOA,F) 

COM-H).'' RJ( 150) .RV( 15 0) 

COMPLEX F ,N ,mH ,0,00 ,P ,IIK ,AI 

pix)=cexp ( CM PLX ( .o.x ) ) • 

PI=3.1A1592 

PH = PHn=i=P t /I RQ. ^ 

B = 2 . =?PI /XLMO A 

XM = 1 ,3 

AI =CriPLX ( .0 , 1. ) 

iM= ( .0, .0) ^ 

DO 5 M=l,20 

I=M-1 ^ 

EM =2. 

1F(1 .80. 0)FM = 1. 

RM = Pe SSJ( ARC, , 1 /XM ,pj ,BY , I6R ) 

MM=EM*RM»P ( I *PI / ( 2 .»XM ) ) »COS( I »PH/XAh) 

M=N+NM 

AMM = CABS(MM ) 

IFtA.XM .LT. l.E-30)0n TO AA 

5 CON T I MU E 

4A 0 = ( . 0 , • 0 ) 

no 13 M=i ,20 

I=M-1 

RM = RE5SJ( ARC, (2 .=1=1 +1. )/ ( 2.»XM 1 ,BJ, BY ,1ER ) 

GM = Rh*P ( ( 2 .-I +1 . ) / ( A ,*XM 1 ) ACOSI ( 2.A1+1. )*PH / ( 2.AXi\i ) ) 

0=G-t-GM 

AGM = CABS(GM ) 

1F(AGM ,LT. l.E-30)G0 TO 55 

15 COMTIMUE 

55 C=SOHT( 2,/(Pl A-R ) )/XM ^ 

G = GAC«P (-P I /A . ) 

M=^»CAP(-P1/A, )/2, 

R0 = R6AL(UK ) 

A1U = AIMAG(UK ) ^ 

IF ((PHI) .6(3. 90.) .OR. (PHI) .60. 27 0. ))GU TO 1 1 

1F((R(I .60. 0.) .AMD. (All) .60. 0.))G0 TO 22 

F = ( AI ACOS(PHn=N+UKAG)/( Ap»COS(PH)+UK) 

on TO 33 • 

11 F=G 

on TO 33 

22 F=N 

33 RETURN 

EMO 
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